Generation of tripartite Einstein–Podolsky–Rosen steering by cascaded nonlinear process
Liu Yu1, Liang Su-Ling1, Jin Guang-Ri1, Yu You-Bin1, †, Lan Jian-Yu2, He Xiao-Bin2, ‡, Guo Kang-Xian3
Key Laboratory of Optical Field Manipulation of Zhejiang Province, Department of Physics, Zhejiang Sci-Tech University, Hangzhou 310018, China
Shanghai Institute of Space Power Sources, Shanghai Academy of Spaceflight Technology, Shanghai 200245, China
School of Physics and Electronic Engineering, Guangzhou University, Guangzhou 510006, China

 

† Corresponding author. E-mail: ybyu@163.com 33349445@qq.com

Project supported by the National Natural Science Foundation of China (Grant Nos. 61975184, 91636108, and 61775043), the Natural Science Foundation of Zhejiang Province, China (Grant No. LY18A040007), the Science Foundation of Zhejiang Sci-Tech University (Grant Nos. 19062151-Y and 18062145-Y), and the Open Foundation of the Key Laboratory of Optical Field Manipulation of Zhejiang Province, China (Grant No. ZJOFM-2019-002).

Abstract

A scheme is proposed to generate genuine tripartite Einstein–Podolsky–Rosen (EPR) steering in cascaded nonlinear process of the fourth-harmonic generation. The second-harmonic is generated by the first double-frequency process in an optical superlattice. Then, the fourth-harmonic is produced by the second cascaded double-frequency process through quasi-phase-matching technique in the same optical superlattice. The genuine tripartite EPR steering among the pump, the second-harmonic, and the fourth-harmonic beams can be obtained by this cascaded nonlinear process according to a criterion for genuine multipartite quantum steering. The quantum steering properties are discussed by adjusting the parameters related to the cascaded nonlinear system. The present research provides a reference scheme and data for obtaining good multipartite EPR steering in experiment and can advance the applications of quantum steering in the quantum information processing.

1. Introduction

The issue of quantum correlation has attracted a lot of interest. The interpretation of quantum correlation can be dated back to Einstein, Podolski, and Rosen’s research of the completeness of quantum mechanics in 1935.[1] This research is the famous Einstein–Podolski–Rosen (EPR) paradox. In response to EPR paradox, Schrödinger proposed the concept of quantum steering (EPR steering)[24] which described by using different observables to detect one of the particles, causing the corresponding one to collapse to a different state. In the decades after this concept was proposed, scientists have been focusing on quantum entanglement[58] and Bell non-locality.[9,10] Until 2007, Wiseman et al.[11,12] found that in the form of quantum information tasks, quantum steering can predict a new property that cannot be described by the local hidden state. Wiseman et al.[11] gave the rigorous definition of quantum steering which lies between quantum entanglement and Bell non-locality. Because of the feature of asymmetry, EPR steering has been deep researched in recent years which aims to a better understanding of the steerability.

The criterion of quantum steering is the theoretical basis for realizing EPR steering, it has obtained the most research results recently.[1319] The simplicity, operability, and practicality of its criteria will greatly facilitate and stimulate the realization and development of EPR steering. In theory, various forms of EPR steering criteria have been developed and applied to EPR steering experiments of different systems and different ideas. In 1992, Ou et al.[20] achieved the violation of EPR steering inequalities for continuous variable by using a nondegenerate optical parametric oscillator (OPO). In 2012, Smith et al.[21] closed the detection of loophole in the optical system and observed the phenomenon of EPR steering. In the same year, Wittmann et al.[22] further closed the free loophole in the experiment. In addition, there are experimentally operable EPR steering criteria.[14,23] For continuous variables, the Reid criterion can be used to detect EPR steering.[15] In 2013, He and Reid put forward a criterion of genuine multipartite EPR steering.[24] Hereafter, by using a continuous variable optical system, a scheme of quantum security communication was proposed and finished.[25] Recently, it has been proved theoretically that EPR steering is an asymmetric quantum nonlocality[26,27] and has been demonstrated in experiments.[28,29] In 2017, Olsen theoretically studied the quantum correlation in the OPO with an injected signal and found the asymmetric steering is controllable.[30] Subsequently, bipartite entanglement and EPR steering correlation was shown by studying the cascaded χ(2) system[31] and the cascaded third-harmonic generation process.[32] Three schemes were experimentally demonstrated to manipulate the direction of EPR steering.[33] EPR steering was proposed can be generated via atomic coherence,[3436] optomechanical systems,[37] and cascaded nonlinear processes in optical superlattice.[38] These great achievements have made us understand deeper about the EPR steering. So far, the EPR steering is fast becoming a key instrument in secure quantum communication[3941] and quantum key distribution.[42] Tripartite entanglement in cascaded nonlinear process of the fourth-harmonic generation was investigated both in an optical cavity[43] and without optical cavity.[44] However, multipartite quantum steering in this cascaded nonlinear process has not been investigated up to now.

In this paper, we present the genuine tripartite quantum steering in cascaded nonlinear process of quasi-phase-matching fourth-harmonic generation in an optical cavity according to the criterion of genuine multipartite EPR steering which put forward by He and Reid.[24] We also discuss the quantum steering properties by adjusting the parameters related to the system. Our research provides a reference scheme and data for obtaining good multipartite quantum steering in experiment and can advance the applications of quantum steering in quantum information processing. The rest of this paper is arranged as follows. In Section 2, the equations of motion are solved and the stationary solutions are obtained in the positive-P representation.[45] The genuine tripartite EPR steering is presented and the quantum characteristics is discussed in Section 3. Finally, a brief summary about this work is given in Section 4.

2. Equations of motion and the stationary solutions

In this work, A fundamental field with the frequency ω0 is incident onto the optical cavity which an optical superlattice (OSLT) is placed inside which can be seen in Fig. 1(a). The second-harmonic field at frequency ω1 is generated by the cascaded nonlinear sum-frequency process, where ω1 = 2ω0. Then, by a second cascaded process between the second-harmonic and second-harmonic fields, the fourth-harmonic field at frequency ω2 = 2ω1 = 4ω0 is generated. In the two cascaded nonlinear sum-frequency processes, the phase mismatching can be simultaneously compensated by the reciprocal vectors G1 and G2 provided by the OSLT through quasi-phase-matching technique[46] as k1 = k0 + k0 + G1 and k2 = k1 + k1 +G2, which is shown in Fig. 1(b).

Fig. 1. (a) Sketch of the optical cavity, (b) quasi-phase-matching schematic diagram for the cascaded nonlinear interaction process.

The interaction Hamiltonian for this cascaded nonlinear process can be written as[43]

where κ0 and κ1 represent the nonlinear coupling constants of optical superlattice, (i = 0,1,2) is the bosonic annihilation operator. Generally speaking, it is taken to be real[47] for the purpose of simplicity. The Hamiltonian of the cavity pumping reads as

where ε is the amplitude of pumping that it is also considered to be real. Following the description of Lindblad superoperator, the losses of the three fields can be written as

where is the system density matrix and γi (i = 0,1,2) represent the cavity loss at ωi. And the master equation of this system can be expressed as

One can map the master equation onto Fokker–Planck equation (FPE) in the positive-P representation for the sake of being able to study quantum steering characteristics better. Here, αi and are independent variables which correspond to αi and when the averages of products converge to normally ordered operator expectation values.[30] Thus, we obtain the FPE of the system as

Following the standard processing, the equations of motion can be derived in the positive-P representation,

where ηi(t) (i = 0,1,2) are the complex Gaussian noise terms which are different from the case in Ref. [43]. They satisfy the relations , , and . In these equations, α0, α1, and α2 correspond to the fundamental, second-harmonic, and fourth-harmonic field amplitudes in the positive-P representation, respectively. By setting dαi/dt = 0, we obtain the steady-state solutions when the noise terms are ignored and found the steady-state solution A1 of α1 satisfies the following relationship

where J0 = 2κ0A1 + γ0, , and J2 = ε κ0 γ2. And we take γ1 = γ2 = γ for the purpose of simple calculations. However, equation (7) is difficult to solve analytically since it requires the solution of a fifth-order polynomial. For this reason, the other two steady-state solutions A0 and A2 can be represented with A1 as follows:

In the following, we can expand the variables into their steady-state expectation values and small Gaussian fluctuations terms close to the steady-state values as αi = Ai + δ αi (i = 0,1,2). Based on that, by means of linear processing approach, equation (6) can be rewritten as

We can rewrite it in a matrix form as

with

where B is the matrix of noise terms contains the steady-state solutions which found from the original FPE and dW is a vector of Wiener increments.[45] The drift matrix A can be written in the following form:

Only when the above drift matrix A has no negative eigenvalues, can the system be in a steady state. In Fig. 2, we show that the real parts of eigenvalues of the drift matrix A (RPEA) versus (a) ε, (b) γ0, (c) γ, and (d) κ1/κ0, respectively. As figure 2 shows, the eigenvalues of the drift matrix A are all non-negative in parameter ranges. The satisfaction of this condition, i.e., the drift matrix A has no negative eigenvalues, guarantees that we can discuss quantum steering characteristics in the steady state. The difference in the selected parameters in the present work in order to obtain better quantum steering causes the result to be different from our previous work[43] for obtaining better quantum entanglement state. Then, the satisfaction of the condition of matrix A allows us to calculate the intracavity spectra in the foundation of Fourier transformation via

where ω is the Fourier analysis frequency, I is the identity matrix, and D = B BT is a 6 × 6 matrix. According to the standard input–output relationship,[49] the output spectra can be calculated from Eq. (13).

Fig. 2. Real parts of eigenvalues of drift matrix A (RPEA) versus (a) ε, (b) γ0, (c) γ, and (d) κ1/κ0, respectively.
3. Genuine tripartite EPR steering

There are many criteria for EPR steering, such as the criteria of Reid[15] in continuous variable, the criteria of multipartite system,[24] and experimentally operable EPR steering criteria.[14,50] Bipartite asymmetric quantum steering in the cascaded χ(2) system[30,31] and the third-harmonic quantum steering[32] are studied, there is almost no genuine EPR steering involved. In this paper, the properties of quantum steering among three fields in the cascaded sum-frequency process will be investigated according to the criterion of multipartite system.[24] We give the quadrature definitions as αi = Xi + iYi, and , all the relationships we present involved Xi and Yi. Thus, based on the criterion of multipartite system, we write the three equalities in the form as

By measuring the values of Si (i = 0,1,2) and comparing with 1 to confirm whether the system i is steerable or not. When the condition Si < 1, it means the success of EPR steering of system i and also indicates that the three fields are bipartite quantum steering each other. There is no doubts that when the measurement of

will be sufficient to demonstrate genuine tripartite steering.[24]

Parameters will affect the results that whether the system is steerable or not in the cascaded sum-frequency process, thence we take the normalized analysis frequency Ω = ω/γ0 into account first. The values of Si and Stot versus the normalized analysis frequency Ω with γ0 = 0.001, γ1 = γ2 = 0.03, κ0 = 0.05, κ1 = 1.5κ0, and ε = 0.06κ0, respectively, is shown in Fig. 3. The choice of parameters is partially identical to our previous work[43] in order to obtain better tripartite quantum steering. As shown in Fig. 3, Si (i = 0,1,2) are all below 1 in a wide range of the analysis frequency which indicates that bipartite quantum steering exists in the three fields. And it can also be found that the value of Stot < 1 across the whole range of analysis frequency Ω from Fig. 3. That is to say, the genuine tripartite steering condition is satisfied, which further confirms the genuine tripartite steering.

Fig. 3. The values of S0 (red dashed line), S1 (black dot-dashed line), S2 (blue thin line), and Stot (green thick line) versus / the the normalized analysis frequency Ω with γ0 = 0.001, γ1 = γ2 = 0.03, κ0 = 0.05, κ1 = 1.5κ0, and ε = 0.06κ0, respectively.

In quasi-phase-matching fourth-harmonic generation process, we can adjust nonlinear coupling parameters. Figure 4 plots the values of Si and Stot versus the nonlinear coupling parameter κ1/κ0 with γ0 = 0.001, κ0 = 0.03, γ1 = γ2 = 0.05, and ω = 4γ0, respectively. We can clearly see that the values of the Si and the Stot increase as the increase of the ratio κ1/κ0 in this range. The value of the Stot will exceed 1 when the ratio of κ1/κ0 is larger than 12 which proves that it is not a genuine tripartite steering. Only when the condition Stot < 1, will the genuine tripartite steering be demonstrated. We found that the value of the nonlinear coupling parameter κ0 directly affects the value of Si and Stot which determined whether the system is steerable and genuine tripartite steering or not in the cascaded sum-frequency processes.

Fig. 4. The values of S0 (red dashed line), S1 (black dot-dashed line), S2 (blue thin line), and Stot (green thick line) versus / the nonlinear coupling parameter κ1/κ0 with γ0 = 0.001, κ0 = 0.03, γ1 = γ2 = 0.05, and ω = 4γ0, respectively.

Figure 5 depicts the values of Si and Stot versus the damping rates γ1/γ0 with γ0 = 0.001, κ0 = 0.01, κ1 = 2κ0, ε = 0.03κ0, and ω = 5γ0, respectively. One can see that the value of Si and Stot are far less than 1 comprehensively. It is below 1, although the values of Si and Stot reache the peak at a rate of about 1.5, which indicates that bipartite EPR steering and genuine tripartite steering can be generated in the cascaded nonlinear process.

Fig. 5. The values of S0 (red dashed line), S1 (black dot-dashed line), S2 (blue thin line), and Stot (green thick line) versus / the damping rate of γ1/γ0 with γ0 = 0.001, κ0 = 0.01, κ1 = 2κ0, ε = 0.03κ0, and ω = 5 γ0, respectively.
4. Conclusion

In this paper, the cascaded nonlinear process of the fourth-harmonic generation is analyzed in the positive-P representation. We demonstrate that the genuine tripartite quantum steering can be obtained in the cascaded nonlinear process based on the criterion for the genuine multipartite quantum steering.[24] The parameters that affect the results are also discussed and we found the better quantum steering can be obtained when the cavity loss is lower. The value of the nonlinear coupling parameter directly determined whether the system is steerable or not. In general, quantum steering, as an emerging quantum information processing sub-direction, has some important and valuable results, but there are still many unknowns waiting for further study. The tripartite steering can also be generated by four-wave mixing (FWM) in rubidium atoms combined with a linear beamsplitter or cascaded a second FWM and their results, showing that the cascaded FWM scheme is a promising candidate to generate stronger multipartite EPR steering.[37] Compared with their scheme, our scheme only needs one pump and one optical superlattice, which is more simple in experiment. However, their scheme has more flexibility to manipulate the monogamy relation. In addition, in their scheme, the correlation level of three-mode is weaker than the case of twin beams since the three-mode correlation condition is more stringent,[51,52] which is consistent with our results as shown in Figs. 25. In previous single-pass cascaded nonlinear experiments, the frequency conversion efficiencies are up to 20% for the second-harmonic generation and close to 23% for the cascaded third-harmonic generation by quasi-phase-matching technique.[46,53,54] Theoretically, the frequency conversion efficiency for the cascaded fourth-harmonic generation is up to 4%. In our scheme, when taking it into optical cavity, the conversion efficiency will be further increased. Therefore, we think that our scheme is feasible in experiment. The difficulty is how to design an optical resonator so that three light fields of different frequencies can resonate at the same time. We hope the results reported in this paper could enable researchers to understand the EPR steering deeper, and give some help for potential applications in quantum information processing.

Reference
[1] Einstein A Podolsky B Rosen N 1935 Phys. Rev. 47 777
[2] Schrödinger E 1936 Math. Proc. Cambridge Philos. Soc. 32 446
[3] Schrödinger E 1935 Math. Proc. Cambridge Philos. Soc. 31 555
[4] Li H Z Han R S Zhang Y Q Chen L et al. 2018 Chin. Phys. 27 120304
[5] Liu Y H Wu L Yan Z H Jia X J Peng K C et al. 2019 Acta Phys. Sin. 68 034202
[6] Ma Y Y Feng J X Wan Z J Gao Y H Zhang K S et al. 2017 Acta Phys. Sin. 66 244205
[7] Zhao C Y Zhang C M 2018 Chin. Phys. 27 084204
[8] Yu Y B Wang H J Feng J X 2011 Chin. Phys. Lett. 28 090304
[9] Bell J S 1964 Physics (Long Island City, N.Y.) 1 195
[10] Liao C G Chen Z H Luo C L 2010 Acta Phys. Sin. 59 8526
[11] Wiseman H M Jones S J Doherty A C 2007 Phys. Rev. Lett. 98 140402
[12] Jones S J Wiseman H M Doherty A C 2007 Phys. Rev. 76 052116
[13] Bennet A J Evans D A Saunders D J Branciard C Cavalcanti E G Wiseman H M Pryde G J et al. 2012 Phys. Rev. 2 031003
[14] Cavalcanti E G Jones S J Wiseman H M Reid M D 2009 Phys. Rev. 80 032112
[15] Reid M D 1989 Phys. Rev. 40 913
[16] Evans D A Cavalcanti E G Wiseman H M 2013 Phys. Rev. 88 022106
[17] Evans D A Wiseman H M 2014 Phys. Rev. 90 012114
[18] Skrzypczyk P Cavalcanti D 2015 Phys. Rev. 92 022354
[19] Zhen Y Z Zheng Y L Cao W F Li L Chen Z B Liu N L Chen K et al. 2016 Phys. Rev. 93 012108
[20] Ou Z Y Pereira S F Kimble H J Peng K C 1992 Phys. Rev. Lett. 68 3663
[21] Smith D H Gillett G de Almeida M P et al. 2012 Nat. Commun. 3 625
[22] Wittmann B Ramelow S Steinlechner F et al. 2012 New J. Phys. 14 053030
[23] Saunders D J Jones S J Wiseman H M Pryde G J 2010 Nat. Phys. 6 845
[24] He Q Y Reid M D 2013 Phys. Rev. Lett. 111 250403
[25] Armstrong S Wang M Teh R Y et al. 2015 Nat. Phys. 11 167
[26] Bowles J Vétesi T Quintino M T Brunner N 2014 Phys. Rev. Lett. 112 200402
[27] Quintino M T Vértesi T Cavalcanti D et al. 2015 Phys. Rev. 92 032107
[28] Hädchen V Eberle T Steinlechner S et al. 2012 Nat. Photon. 6 596
[29] Wollmann S Walk N Bennet A J Wiseman H M Pryde G J 2016 Phys. Rev. Lett. 116 160403
[30] Olsen M K 2017 Phys. Rev. Lett. 119 160501
[31] Olsen M K 2017 Phys. Rev. 96 063839
[32] Olsen M K 2018 Phys. Rev. 97 033820
[33] Qin Z Z Deng X W Tian C X et al. 2017 Phys. Rev. 95 052114
[34] Zhong W X Cheng G L Hu X M 2017 Opt. Express 25 11584
[35] Wang L Lv S C Jing J T 2017 Opt. Express 25 17457
[36] Liu Y Cai Y Xiang Y Li F Zhang Y P He Q Y et al. 2019 Opt. Express 27 33070
[37] Li J Zhu S Y 2017 Phys. Rev. 96 062115
[38] Liu Y Liang S L Jin G R Yu Y B 2020 Opt. Express 28 2722
[39] He Q Y Rosales-Zrate Adesso L G Reid M D 2015 Phys. Rev. Lett. 115 180502
[40] Xiang Y Kogias I Adesso G He Q Y 2017 Phys. Rev. 95 010101
[41] Chiu C Y et al. 2016 njp Quantum Inf. 2 16020
[42] Walk N Hosseini S Geng J et al. 2016 Optica 3 634
[43] Yu Y B Wang H J 2011 J. Opt. Soc. Am. 28 1899
[44] Yu Y B Wang H J Zhao J W Ji F M Wang Y J et al. 2016 Laser Phys. Lett. 13 085203
[45] Gardiner C W 1991 Quantum Noise New York Springer 203
[46] Zhu S N Zhu Y Y Ming N B 1997 Science 278 843
[47] Ferraro A Paris M G A Bondani M Allevi A Puddu E Andreoni A et al. 2004 J. Opt. Soc. Am. 21 1241
[48] Yu Y B Wang H J Xiao M Zhu S N 2011 Opt. Express 19 13949
[49] Collett M J Gardiner C W 1984 Phys. Rev. 30 1386
[50] Saunders D J Jones S J Wiseman H M Pryde G. J 2010 Nat. Phys. 6 845
[51] Li C B Jiang Z H Zhang Y Q Zhang Z Y Wen F Chen X X Zhang Y P Xiao M et al. 2017 Phys. Rev. Appl. 7 014023
[52] Zhang D Li C B Zhang Z Y Zhang Y Q Zhang Y P Xiao M et al. 2017 Phys. Rev. 96 043847
[53] Zhu S N Zhu Y Y Qin Y Q Wang H F Ge C Z Ming N B et al. 1997 Phys. Rev. Lett. 78 2752
[54] Qin Y Q et al. 1998 Appl. Phys. Lett. 84 6911